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Abstract We undertake a hydrodynamical study of a mixture of tightly cou¬ 
pled primordial radiation, neutrinos, baryons, electrons and positrons, together 
with a gas of already decoupled dark matter WIMPS and an already existing 
“frozen” magnetic field in the infinite conductivity regime. Considering this 
cosmic fluid as the source of a homogeneous but anisotropic Bianchi I model, 
we describe its interaction with the magnetic field by means of suitable equa¬ 
tions of state that are appropriate for the particle species of the mixture be¬ 
tween the end of the leptonic era and the beginning of the radiation-dominated 
epoch. Fulfilment of observational bounds on the magnetic field intensity yields 
a “near FLRW” (but strictly non-perturbative) evolution of the geometric, 
kinematic and thermodynamical variables. This evolution is roughly compara¬ 
ble to the weak field approximation in linear perturbations on a spatially flat 
FLRW background of sources in which the frozen magnetic fields are coherent 
over very large supra-horizon scales. Our approach and results may provide in¬ 
teresting guidelines in potential situations in which non-perturbative methods 
are required to study the interaction between magnetic fields and the cosmic 
fluid. 
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1 Introduction 

A full account of the origin and evolution of cosmic and astrophysical magnetic 
fields is a long standing problem in contemporary theoretical physics. Magnetic 
fields have been observed in all scales: our own galaxy [1] , low [2] and high red- 
shift [3] galaxies, and up to scales of galaxy clusters [4] and superclusters [5] (for 
comprehensive reviews see 0)- There is also indirect evidence, from gamma- 
ray observations of blazars [7] of coherent intergalactic magnetic fields in low 
density regions. While magnetic fields can result as a consequence of a variety 
of local astrophysical processes involving electric charges in motion (accretion 
into active galactic nuclei or compact objects, ionization of intergalactic gas 
and interactions between baryons and Cosmic Microwave Background (CMB) 
photons), all this multi-scale observational evidence supports the hypothesis 
that (at least) some of them may have a primordial origin from “seed” mag¬ 
netic fields that emerged during some early evolution stage (or stages) of the 
cosmic fluid. As shown extensively in recent comprehensive topic reviews 
[iniin] on cosmological magneto-genesis, seed magnetic fields can be generated: 

— during inflation: by breaking the conformal invariance of Maxwell theory 
through the interaction of electromagnetic fields with dilaton-like scalar 
fields or axion-like pseudoscalar field |12j . or by assuming conformal in¬ 
variance with an FLRW background with negative spatial curvature |13j : 

— in the various phase transitions after inflation by non-standard couplings 
that may generate helical magnetic fields on small scales (of the order of 
the horizon scale at these transitions). These fields can be amplified to 
all scales through an inverse cascade from the conservation of magnetic 
helicity M- Magnetic fields can also emerge through plasma mechanisms 
(Weibel effect) during QCD phase transition [TK] : 

— during the classical plasma phase before recombination |16| . by assuming 
that the photon-electron coupling in Compton scattering is much stronger 
than the coupling between photons and ions, leading to a density fluc¬ 
tuation field generating magnetic fields through pressure anisotropy and 
induced electric currents associated with the differences in the bulk elec¬ 
tron and proton velocities. 

In the above mentioned articles on early times cosmological magneto-genesis 
the dynamics of the magnetic interaction in the cosmic fluid is examined by 
means of a variety of techniques and theoretical frameworks depending on 
the specific transition: numerical simulations. Kinetic Theory, Statistical Me¬ 
chanics, or gauge invariant vector perturbations over a spatially flat FLRW 
background (see details in [Hl IHlUOlllI] , see other approaches to primordial mag¬ 
netism in pTlfT^fT^ b However, a classical Maxwellian electromagnetic interac¬ 
tion (not necessarily restricted to early cosmic times) can also be incorporated 
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into the full dynamics of General Relativity through an elegant first order 
system based on covariant objects defined in a 4-velocity frame (the 1+3 for¬ 
malism [nilinilll])- In particular, this formalism is useful to study (following a 
perturbative or non-perturbative approach) the dynamics of a “frozen” mag¬ 
netic field without electric currents through the limit of infinite conductivity 
Whenever a spatially flat FLRW background is assumed 
in a perturbative treatment, a “weak field approximation” can be defined in 
which the anisotropic stresses of the the magnetic field (associated with vector 
perturbations) are neglected and the latter field is described (at first order) 
as a sort of scalar source that dilutes as R ^ 

While this approximation is not gauge invariant, it leads to the same results 
(see proof in [26) 1 as gauge invariant perturbations when infinite conductivity 
is assumed. This weak field approximation was used also in earlier literature 
p [?fl[^ (see update in |^[5D] 1 in a purely FLRW context to examine the 
magnetic interaction at early cosmic times. 

A different approach to cosmological magnetism that is non-perturbative 
(but restrictive) consists in considering a fluid (or fluid mixture) together with 
a magnetic field (in the infinite conductivity approximation) as sources of 
exact solutions of Einstein’s equations. Work along these lines is found in older 
literature (see pioneering work in [31 j l considering very simple highly idealised 
fluids (even dust) in solutions with cylindrical symmetry |32j , or homogeneous 
but anisotropic Bianchi models |33j and, in particular, spatially flat Bianchi I 
models I51IM1I3S], which provide the simplest spacetime geometry compatible 
with the anisotropic stresses that characterise the magnetic interaction. More 
recently a Bianchi I model was used to examine how the anisotropic 
stresses produced by free streaming neutrinos (described in terms of relativistic 
Kinetic Theory) allow for the fulfilment of the late time constraints placed by 
the CMB on the anisotropy of large scale (supra-horizon) magnetic field^ 

In the present article we examine a magnetised early Universe cosmic fluid 
mixture as a source of a Bianchi I model. Instead of considering idealised 
fluids whose interaction with the magnetic field is unspecified (as in [^[3^ 
[MllSilISS]), or a collision-less kinetic theory approach (suited to free stream¬ 
ing collision-less conditions as in m), we consider fluid sources that satisfy 
physically motivated equations of state that are well suited for an early Uni¬ 
verse cosmic fluid interacting with a magnetic field within a hydrodynamical 
regime, namely: a mixture of ideal gasses of magnetised fermions [39ll40ll4T] 
whose equation of state reflects the full anisotropic effects of the magnetic field. 
Evidently, this field introduces anisotropic momentum fluxes that modify the 
energy-momentum tensor and thus affect the evolution of the state variables. 
Hence, the anisotropic pressure terms associated with the magnetic field nec¬ 
essarily contain a (“pure”) classical Maxwell term [42], but must also modify 
the equations of state of the fluid sources. We remark that in previous work we 

^ We remark that Bianchi I models have received renewed attention m> since present 
observational bounds from Planck and WMAP cannot rule out (in principle) the possibility 
that the large scale structure of the Universe is not perfectly isotropic from a statistical 
point of view. 
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have studied the gravitational collapse of such magnetised sources in a Bianchi 
I geometry, considering the case of zero [4311^1^ and finite temperatures [46] . 

In order to examine the hydrodynamical evolution of the magnetised fluid 
mixture described above, we assume the magnetic field to be the only source of 
anisotropy of the energy-momentum tensor and consider two separate cosmic 
epochs where the effect of the magnetic field can be particularly relevant for 
a Universe full of free electric charges: (i) the end of the leptonic era, and: 
(ii) the beginning of the radiation dominated era immediately before cosmic 
nucleosynthesis. Although all fermions, neutral and charged, interact with the 
magnetic field, we will assume an equation of state in which the contribution 
of electrons and positrons is dominant. This is a good approximation, since in 
the cosmic times we are interested the constraint eB ^ „ holds and the 

interaction of the magnetic field with protons and neutrons can be neglected 
from the statistical properties of these particles. 

Besides introducing anisotropic stresses, primordial magnetic fields lead to 
various important effects on the evolution of early Universe sources. Hence, 
stringent bounds need to be imposed on their field strength in order to comply 
with numerous observational constraints, including limits of the abundances of 
small anisotropic stresses from older literature based on the COBE four-year 
anisotropy data [47ll48j . In particular, we are concerned on the effects of the 
magnetic field on cosmic nucleosynthesis, since it is the cosmological event in 
the cosmic stages under examination that is “closest” to us. Magnetic fields 
lead to various different effects on cosmic nucleosynthesis: (i) their contribution 
to the energy density content of the cosmic fluid affects the expansion rate; 
(ii) the electron-positron quantum statistics is modified, and so is the rate of 
neutron beta decay (see [8| for a review). Considering these effects together, 
an upper bound of (Bq) < 3 x 10“'^G is necessary at length scales of the order 
of the Hubble horizon size at BBN time An updated value of this 

bound is: {Bq) ^ 1.5 x 10“®G (this is related to the local field amplitude B 
contributed from all wavelengths) |29j . These bounds were obtained without 
considering the anisotropy introduced by the magnetic field. On the other 
hand, fitting GMB observations leads to a stronger constraint that appears to 
exclude homogeneous cosmological magnetic fields much stronger than 10“® 
G 0137]. 

Besides considering observational bounds on the magnetic field, we also 
discuss how a comparison of our approach and results can be made with those 
based on cosmological perturbations discussed and summarised in pnniiii] 
and in previous work based on the weak field approximation 
This comparison is facilitated by our choice of a Bianchi I model that be¬ 
comes a spatially flat FLRW model when the magnetic field vanishes {i.e. the 
magnetic field is the sole cause of anisotropy of the source), as in this case a 
sufficiently weak field can be regarded approximately as a sort of supra-horizon 
magnetic perturbation on a spatially flat FLRW background (see comments 
on this issue in 137 ]), but it lends naturally with a comparison in the weak 
field regime studied in [2511261130] and in |8]. Without neglecting the impor¬ 
tance of observational bounds, the obtained results may be useful in potential 
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situations in which strong cosmic magnetic field need to be examined in a 
non-perturbative manner. 

The paper is organised as follows: in Sec. we examine the Einstein- 
Maxwell field equations that govern the dynamics of a magnetised cosmological 
fluid mixture in a Bianchi I geometry. The fluid sources and equations of state 
for the constituents of this magnetised fluid are discussed in detail in Sec. 
The Einstein-Maxwell system is transformed into a first order system given in 
terms of dimensionless variables in Sec. 4. The perfect fluid and FLRW limits 
are examined in Sec. 5, while in Sec. 6 we discuss the possible comparison with 
linear perturbations and the weak field regime and the late time evolution 
of the model. In Sec. we report our numerical results for the evolution of 
relevant thermodynamical and dynamical parameters. Finally, we provide our 
conclusions in Sec. HI 


2 Einstein Maxwell equations for a Bianchi I model 

Considering rectangular comoving coordinates {t, x, y, z) with x, y, z aligned 
along the principal axes of the shear tensor, Bianchi I models are described by 
the Kasner metricQ 


ds — —dt T OL-^ (t) dx T 0,2 (t) dy -t- cl^ dz . (1) 

In order to describe a matter-energy source that includes magnetic interaction 
for this metric we consider the following energy-momentum tensor 




( 2 ) 


where p = p{t), p = p{t) and are, respectively, the energy density, 

isotropic pressure and the traceless spacelike symmetric part of the stress ten¬ 
sor (a nonzero is necessary to describe a magnetic field, see equations (16) 
([T7|) further ahead) and A is the cosmological constant (whose value is sub- 


and 

jected to observational constraints). The Bianchi I geometry yields a uniform 
magnetic field (constant at each fixed t) of infinite extension along a unique 
spacelike tetrad vector direction that can always be identified with a space 
axis (say, the z-axis). Hence, in a comoving frame with 4-velocity this 

stress-energy tensor takes the diagonal form: 


n = diag [-P - A, P^- A, P^ - A, P|| - H] , (3) 

where P± = P±(t) and Py = P\\it) are the pressures {i.e. eigenvalues of the 
tensor PI^v) in the directions perpendicular and parallel to the magnetic field 
in the frame defined by the coordinates (t, x, y, z). 

The dynamics of a cosmological fluid allowing for a magnetic interaction 
follows from the coupled Einstein-Maxwell field equations (G^j/ = SttT^,^), 


2 


Unless specified otherwise, we use natural units G = c = 1. 
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which for the Bianchi I geometry Q take the form: 


G X _ 


-G^ = 


0203 

0203 

0103 

O1O3 

0102 


21 + 51 = -8,(Px-/1), 

02 03 

Ol O3 


_G/ = + ^ + ^ = -87r(P|| - 2 I), 

0102 Ol O2 
0102 0103 0203 


-G/ = 


0102 0103 


+ =87r(p + 2l). 


02 O 3 


(4) 

(5) 

( 6 ) 
(7) 


where d = u°‘a^a = o,t. The conservation equation -^ = 0 yields: 


P = -(^ + ^)(P3_+p)-^(P||+p), (8) 

\Ol 02/ 03 


together with the only non~trivial of Maxwell’s equations ;j/ = 0 and 
~ 0) ■ 


Ol 

, 02 

B 


H- 

+ — 

Ol 

02 

B 


(9) 


Notice that B is not a scalar. It can be characterised as the single nonzero 
component of the covariant magnetic held vector in the comoving frame: = 

, where 6^ and are the Kronecker delta and Levi- 

Civita tensors. 


3 Physically motivated field sources 

We will assume as the field sources of the Bianchi I model (i) a tightly coupled 
gas mixture of leptons, baryons and photons, (ii) an already decoupled non¬ 
interacting gas of cold dark matter WIMPS whose contribution is basically the 
rest-mass energy density (whose abundance value can be assumed to be very 
subdominant: 17 cdm 10“®), (iii) an already existing magnetic field B. We will 
follow a purely phenomenological approach to this magnetic field, assuming 
that it emerged from previous fundamental processes whose details are beyond 
the scope of this paper (see brief description on cosmic magneto-genesis in the 
introduction and a comprehensive discussion in pununuEDg 

We will also assume that the magnetic interaction is only significant with 
some of the particle species in (i). The total energy density and pressures for 


^ Assuming a present day value Qq ~ 0.7 for a ACDM background, the contribution of 
the cosmological constant is absolutely negligible in the cosmic times under consideration, 
hence we omit it until we discuss the late evolution of the models in Sec. |6.3[ l 
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the source ©-(i have the form: 


k 

( 10 ) 

Pa_=Px +P^ + P'' + ^Pt, 

k 

( 11 ) 

k 

( 12 ) 


where k runs over all particles that could interact with the magnetic field, and 
the upper indices B, 7 and respectively denote the density and pressures 
associated with the magnetic field, the photons and the neutrinos. 

Because of the symmetries of the Bianchi I model the magnetic field in a 
comoving frame is necessarily homogeneous {i.e. purely time dependent) and 
(as we have chosen it) to be everywhere oriented along the z-direction. The 
magnetic density and pressures are then 

Pf = -Plf = = p^, (13) 

leading to a positive pressure term along the x and y-axes and a negative 
pressure along the field direction. This negative pressure can be interpreted as 
a tension or elasticity of the field lines, which tend to remain as “straight” as 
possible by reacting to any effect that distorts them [4^150] . Consequently, an 
anisotropic pressure must give rise to an anisotropic expansion law. 

Photons contribute to the source with an isotropic pressure and energy 
density given by the equation of state 

p^T) = pyr)/3, p^{T) = a^T\ (14) 

where Og = 7r^/15 and T is the temperature. On the other hand the isotropic 
pressure and energy density of neutral leptons (neutrinos and antineutrinos) 
are given by [M] 


p^T) = p^T)/3, p''(T) = 7yagTVl6, (15) 

with g' equals to 2 x 3, hence we have taken into account the three different 
species of neutrinos and antineutrinos. 

It is safe to neglect the interaction of photons, neutrinos and WIMP’s with 
the magnetic field. However, leptons and baryons could interact with the field 
through their charges (if they are charged) and via their anomalous magnetic 
moment (if they are neutral). In any case the interactions with the field (via 
charges or anomalous magnetic moment) lead to a momentum-energy tensor 
with anisotropic stresses We assume hereafter that the only particle 

species interacting with the magnetic field are electrons and positrons (charged 
leptons), protons (charged baryons) and neutrons (neutral baryons). 
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In general the equation of state for these species in presence of a time- 
dependent magnetic field can be written as follows 


P>1 = -[2^ - BM'^, (16) 

= -f2^ (17) 

+TS'^ + fikNK (18) 


where the upper index k denotes generically the electron, positron, proton 
and neutron, = — (dfi^/dB) is the magnetisation, /dT^ 

is the entropy, = — (df2^/d^k) is the particle number density (with /Xfc 
the chemical potential), and 17^ the thermodynamical potential which has two 
contributions: 

= + (19) 

where is the statistical Quantum Field Theory contribution and 

(which does not depend on the temperature and the chemical potential) is the 
well-known Quantum Field Theory vacuum term given by|^ 




1 

47r^ 


OO 

X! / dp\\d^p^ek, 

_i 1 ^ 




( 20 ) 


and the statistical term: 


OO 

' ’ —OO 

( 21 ) 


In the previous equations /3 = l/T, 77 = 1,-1 correspond to the two orienta¬ 
tions of the magnetic moment (parallel and antiparallel) with respect to the 
field, while Sk is the spectrum of the fermions given by: 


Efe 


^Jpl + ‘^\eB\l + ml 




Charged fermions. 
Neutral fermions. 


( 22 ) 


In the equations above rrik denotes the fermion mass and q is the anomalous 
magnetic moment. For charged fermions we need to carry on the following 
substitution: 


/ 


(PpA_ 

(27r)2 


27r 


^(2 - Sio), 
1=0 


(23) 


See definitions in m- The term has non-field-dependent ultraviolet divergencies, 

after renormalisation the Schwinger expression is obtained m- 
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where d{l) = 2 — Sio is the spin degeneracy of Landau levels with I ^ 0. For 
electron-positron pairs we assume a negligible chemical potential {fj, = 0). 
Hence, the last term on the right-hand side of equation (181 vanishes. 


For the temperature of the full duration of the stage of cosmic evolution 
we are interested the protons and neutrons satisfy {eB ^ so 

they contribute to the evolution mostly through their rest energy. Therefore, 
we assume that the only particles that contribute to the anisotropy in the 
pressures are the electrons and positrons. 


4 Einstein—Maxwell equations as a first order system 

For a hydrodynamical numeric framework it is necessary to transform the 
second order Einstein-Maxwell equations Q-Q into a first order system of 
evolution equations for local kinematic covariant objects (see Appendix 0 
and the relevant state variables. We also need to introduce the following di¬ 
mensionless evolution parameters: t, H, S 2 and S 3 given by 


d Id 
dr Hit dt ’ 

e/3 


n = 


Hit 


S 2 


T = Hit(t- ti), 

_ ^2 C _ ^3 
“ rr ’ *^3 — „ j 

-tl ^ -tl ^ 


(24) 

(25) 


where is an inverse length defined by the condition 3H^ = SttGX/c^, where 
A = vriel (with Ac the Compton wavelength and rtif. the electron mass), 

e = -p, is the expansion scalar and E 2 , H 3 are the eigenvalues of the shear 

tensor (see Appendix A). Notice that iJ, roughly corresponds to the Hubble 
length in the outset of the radiation epoch (just before nucleosynthesis). 

We also define for generic magnetised gas mixtures the following dimen¬ 
sionless variables 


II 

II 

—\ 

V ^11 

^11 = 

(26) 


B 

T 

(27) 



r= —, 

rUe 


where p, P±, Py are the total energy density and anisotropic pressures in (10 )- 
(12) and = rn^/e = 4.41xlO^^Gis the critical magnetic field for electrons^ 
Considering the variables (24)-(27), the Einstein-Maxwell field equations be- 


® The critical magnetic field for an electron as defined above is the strength at which 
electron cyclotron energy equals its rest energy. 
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come the following first order system: 


1 


= - ^ (Sf + 2Vj_ +Vii)-l [(^2 + Ssf - S 2 S 3 

5 2 , r =V^- 1P|| - 3HS2, 

53, r = 2 (1P|| - ip_l) - mSs, 

B,r=B{S 3 - 2 n), 

1 


(28) 

(29) 

(30) 

(31) 


T,r = -^ {i2H - S3) [S + V±) + {H + S3) {£ + V\\) - £,bB (£3 - 2H)) , 

(32) 


together with the Hamiltonian constraint 


i£ = -Si - Si - S 2 S 3 + 3■H^ 


(33) 


where £_j- = dS/dT and £ 3 = dSjdB in (32). 

Since we have provided in -([^ the equations of state for all the 
constituents of the total energy density and pressures in (IIo|)-(12), we can 


now provide constraints that link £, H, P||, 'P\_ (or p, B, Py, P±), the differ¬ 
ential equations (28)-(32| and the constraint (33) become a complete and 


self-consistent system whose numerical integration allows us to examine the 
dynamical evolution of a magnetised Universe in the cosmic stages we are 
concerned. Since all thermodynamical functions depend only on B and T (di¬ 
mensionless magnetic field and temperature, respectively), from the numerical 
solutions for these two variables we can obtain the thermodynamical func¬ 
tions. On the other hand, the solutions for S 2 , S 3 and H provide the necessary 
information to study the kinematical evolution of the cosmic fluid. The local 
proper volume can be expressed in terms of B as follows: 


V{t) = U(Ti)exp 3 / Hdr 


and the expression for metric coefficients reads: 


aa(r) = aa(ri)exp (S'a + 'H)dr^, a = 1 , 2 ,3 


(34) 


(35) 


with Si = —{S 2 + S 3 ) (since the shear tensor is trace-free). 


5 The FLRW and perfect fluid limits 

If the magnetic field vanishes {i.e. B —>■ 0) we have P|| = P± and thus 
-t 0 (from ([^, ^ and ([^-([^). As a consequence, ([^ reduces to 
the momentum-energy tensor of a perfect fluid with isotropic pressure: 

— p UflUy -[- P -\- yi Pfiu 5 


P = P||=P^. 


(36) 
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However, an isotropic pressure does not imply an isotropic geometry (i.e. an 
FLRW geometry), as the energy-momentum tensor (361 can still be compatible 
with the inherent anisotropy of the Bianchi I geometry that is present in the 
different fluid expansion rates in (501 and the nonzero shear tensor in (51) 
and (52) (see Appendix [A|) . Hence, a zero magnetic field in a Bianchi I model 
only leads to an FLRW geometry if besides the magnetic field the shear tensor 
vanishes as well. 

The conditions for an evolution with zero or nonzero shear follow readily 
from the evolution equations (29)-(30), which if the source is a perfect fluid 
{V_L = V\\ from (26)) can be integrated formally as 


S'i(t) = S'i(Ti)exp ^-3 JHcIt^ , S2{t) = S'2(Ti)exp ^-3 JHdr^ , 

(37) 

with S 3 = —{Si + S 2 )- Hence, when B = 0 and the energy-momentum tensor 
becomes (36) we can identify the following two possibilities: 


— Evolution with nonzero shear. If the initial values (S'a(Ti) for a = 1,2,3) 
of (at least) one of the shear eigenvalues is nonzero the shear tensor is 
nonzero for all r > t^. Notice that (for example) S 2 {Ti) = 0 implies S' 3 (Ti) = 
—Si{Ti), which is nonzero in general. 

— Evolution with zero shear: if any two of the initial values 5'a(Ti) (for a = 
1,2,3) vanishes, then Si = S 2 = S 3 = 0 holds for all r > r^. The shear 
tensor vanishes. 


We examine below these two limit cases for a zero magnetic field and the 
possibility of studying the latter in a purely FLRW context. 


5.1 Perfect fluid Bianchi I limit 


If H = 0 but the shear tensor is nonzero (Sain) yf 0 for at least one of the S'a) 
we have a Bianchi I model whose source is the same particle mixture described 
in sectionj^but without the magnetic interaction. This source is characterised 
by the equations of state (13), (16)-(18l with B = 0, leading for each particle 
species to 


=p^l= p* = 


pfc =-12'=+TS' + ^iV, (38) 

Pi = -P\\ = = 0, (39) 


where now the statistical thermodynamical potential 1?'= = 12gQp.j,(/i, T) in (21) 


does not depend on B. The remaining equation (20)-(21) remain valid with 


B = 0. The dynamics follows from Einstein’s equations, which is the system 
(28)-(32) and (33) with V± = V\\ = V and without the Maxwell part: i.e. 


wit 


rout (31) since B = 0. 
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5.2 FLRW limit model 

If iJ = 0 and the shear tensor vanishes (S'a(Ti) = 0 for at least two of then 
P± = P\\ = P and S'a = 0 holds for all a = 1 ,2, 3, hence Oi = 02 = 03 = a. It 
is straightforward to show that reduce to the well known spatially flat 

FLRW equations: 


2 d a? K 

2 l^P-) 2 

a 6 


p = -3{p + P)-, 
a 


(40) 


or, in dimensionless form for the particle mixture: 


n, = --(s + r). 


Tr = - 


m {£ + V) 


'.r 


= £, 


(41) 


where V = P/A. 


6 Perturbations, the “weak field” approximation and late time 
evolution 

Magnetic fields are evidently incompatible with a non-perturbed FLRW Uni¬ 
verse, but the expected near homogeneity and high electric conductivity of 
the early Universe cosmic fluid before recombination [iinEniiiiiiii] suggests 
introducing the magnetic interaction through suitable scalar and vector per¬ 
turbations on an FLRW background. Such a perturbative approach of a frozen 
magnetic field is fully justified if the latter is tangled on scales smaller than 
the Hubble horizon m and has lead to a comprehensive literature nnunnn 
in which the perturbations are covariant and gauge invariant. 


6.1 The magnetised Bianchi I model as an “exact” perturbation 


The usage of an exact Bianchi I model to describe a magnetised fluid mixture 
is not, strictly speaking, a perturbative treatment of the magnetic interac¬ 
tion. However, if we assume the anisotropy of the energy-momentum tensor 
to be only caused by a magnetic interaction associated with very weak mag¬ 
netic fields (as required by observational bounds), then the Bianchi I model 
can roughly approximate a “near FLRW” perturbative-like regime that also 
complies with exact spatial flatness (which to great extent holds in the early 
Universe cosmic fluid). Considering this regime is well justified, given the fact 
that deviations from isotropy (in both the matter-energy source and the space- 
time geometry) are expected to be negligible in early cosmic times. Evidently, 
(from @-@, and @) a weak magnetic field implies a near 

perfect fluid: 

\2B^-87rj:, 


SttA 


^ ^ <Cl, 




(42) 
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which implies in turn (from and ([^-([^) approximately equal scale 

factors close to a unique FLRW scale factor: oi « 02 « 03 « a{t), lead¬ 
ing to negligible shear \Sg^\/H = |S'a| ^ 1 and a near isotropic expansion 
y, « a^r/a + 0{Sa), as well as (from ( [M] )) to an approximate “weak field” 
scaling B_r ~ —2B'H + 0{Sa) ^ B ~ a“^. In fact, it is straightforward to 
show that under these conditions the evolution equations (28)-(33l become 
linearised, taking the form of a set of FLRW “background” equations (41) 
plus two extra “first order” linear equations for B and the shear eigenvalues as 
“exact” perturbations (this linearisation of Bianchi I geometry for magnetised 
sources is discussed in [13]). 

However, while this type of “exact” magnetic perturbations can be cast 
in terms of the variables of the standard gauge invariant or the covariant 
formalisms, they are extremely restrictive because the homogeneity of the 
Bianchi model introduces a magnetic field of infinite extent (that permeates 
the whole 3-dimensional space) and it only allows to examine the (small) 
deviation from the FLRW background of the time dependent amplitude of 
the (exactly) perturbed magnetic source, but (lacking spatial dependence) the 
“near FLRW” Bianchi model cannot describe its dependence on spatial scales 
that would arise naturally in the standard perturbation formalisms. In other 
words: a magnetised Bianchi I model that becomes arbitrarily close to its 
FLRW particular case (by considering sufficiently weak magnetic fields) may 
only describe approximately an exact perturbation of a unique scale of infinite 
spatial extent, and as such can approximate only very large scale supra-horizon 
perturbations as, for example, in various string inspired inflationary magneto¬ 
genesis scenarios (see details in [31I31 [TI71[TT] and also in [33]), and in the Bianchi 
I model used in m to study for the interaction of free streaming neutrinos 
and such large scale weak magnetic field. 


6.2 Magnetic fields in an FLRW context 

It is more natural to compare the magnetised Bianchi I model that we have 
examined with previous work in [31l371[331l331[3n] . where (besides considering 
a similar theoretical framework based on Satistical Mechanics) the magnetic 
interaction in the infinite conductivity regime and the cosmic fluid are both 
studied in a purely FLRW context, in which the frozen magnetic field is intro¬ 
duced as part of the background (the “weakly magnetised FLRW” spacetime 
in [30] 1 as a relativistic scalar correction to the energy density and the isotropic 
pressure of a perfect fluid particle mixture considered as source of an FLRW 
metric (see equations (3.26)-(3.29) of [ 8 ] and section IX of [30]). This approach 
represents a further simplification of the weak field limit that is applicable to 
an early radiation dominated era (well before structure formation) in which 
matter and radiation perturbations are also neglected, and can be justified if 
the magnetic held is “not too tangled on scales smaller than the magnetic dis¬ 
sipation scale” (see discussion in [53] and m- This simplihcation of the weak 
held regime is not gauge invariant, but for the cosmic times under considera- 
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tion it yields the same solutions as the weak field (and thus as gauge invariant 
perturbations) in the infinite conductivity regime |25LI26] . 

While the energy density of the perfect fluid mixture described above 


fully coincides with our total energy density in equation ( 10 ), we consider 
anisotropic pressures in order to be consistent with the full non-perturbative 
treatment based on Einstein-Maxwell field equations. Under the approach of 
[siiiaisaiin] the conservation of the frozen magnetic flux yields: 


B 

B 


2 d 


= 0 


B{t) = 


Bit,) 


(43) 


which can be compared with the exact equation ([^ of the Einstein-Maxwell 
system obtained in a Bianchi I Universe (or equation (31) in dimensionless 
variables): 


Ol 02 

Oi 02 


B 

B 


= 0 ^ B{t) = 


B{t,) 

ai{t) a2{t)' 


(44) 


While (43) can be regarded as an approximation of (|44| if anisotropy is neg¬ 
ligible (so that Oi « 02 ), our approach is different from that of [sumuHiiin]: 


we regard the Maxwell equation (44) as part of the coupled Einstein-Maxwell 


system in the Bianchi I geometry, whereas the flux conservation (43) in [51 
r!Tir!5ii?)| merely provides a subsidiary condition for Einstein’s equations in 
the FLRW metric and B is taken as a scalar (which is not correct: it is the 
component of a frame-dependent magnetic field vector). Another important 
difference with respect to liiiTiiiHiiin] is that these authors simply assume 
the FLRW radiation temperature for the magnetised mixture, whereas in our 
Einstein-Maxwell system this temperature needs to be obtained from the evo¬ 
lution equation 


T = - 

P,T 


' Oi 02 / 03 ^ " V Oi 02 ' 


O3 


Ol 


02 


(45) 


which is coupled with the magnetic field and expressed in terms of dimen¬ 


sionless variables (26|-(27) becomes equation (32). However, if we assume an 


FLRW geometry 0 = 01 = 02 = 03 , a perfect fluid source P± = P|| and the ra¬ 
diation equation of state, then ( [4^ becomes the FLRW evolution equation for 
the radiation temperature: T = —Ta/a. Nevertheless, for very weak magnetic 
fields (32) should approximate this radiation temperature evolution law. 

Therefore, the results of [51[?7ll551[55] can always be obtained from our 
Bianchi I based results if we assume the same approximations considered by 
these authors. While these approximations may be well justified for weak fields 
in an early time cosmic mixture, our non-perturbative approach allows us 
to examine magnetised cosmic fluids also when such assumptions cannot be 
justified, as we consider the coupled Einstein-Maxwell system that takes into 
full account the spacetime anisotropy, its effects on the energy density and 
pressures and the vectorial (frame-dependent) nature of the magnetic field. 
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6.3 Late time evolution 

For cosmic times beyond nucleosynthesis we need to consider the evolution of 
the tightly coupled fluid mixture in the continuing radiative classical plasma 
phase before recombination, and then the matter dominated phase after the 
decoupling of radiation and matter, where a hydrodynamical regime is no 
longer a valid assumption. For the parameters and initial conditions we have 
considered (see Sec. the near FLRW evolution at late times yields baryon 
and CDM densities diluting as ~ a“^, radiation, neutrino and magnetic field 
densities as ^ a~^, while the cosmological constant dominates the dynamics 
as the Bianchi I model evolves asymptotically into a ylCDM model. 

Whether we account or neglect possible magneto-genesis or some form 
of magnetic interaction (as discussed in US], see also lanniE]) in the lapse 
between nucleosynthesis and radiation-matter decoupling, the evolution of 
the components of the source we have examined in the subsequent matter 
dominated era needs to be examined by kinetic theory, with the contribution 
of baryons and CDM dominated by their rest mass density and relativistic 
particle species (photons and neutrinos) undergoing a free streaming evolution. 
These are precisely the conditions assumed in m for the a post recombination 
evolution of a Bianchi I model that fulfils present day observational bounds 
on the CMB and tends asymptotically to a ylCDM model. By imposing the 
constraints on the magnetic field that yield the same CMB bounds (see Sec. 
7.4), the Bianchi I model we have examined effectively exhibits the same late 
time evolution as the Bianchi I model of m- In fact, our model also complies 
with the nucleosynthesis constraints on the abundance of light elements (see 



7 Numerical results 

7.1 Stages of cosmic evolution 

We consider for the numerical study of the dynamics of a magnetised Universe 
the following two periods of the cosmic evolution before the nucleosynthesis: 

— Epoch I: End of the leptonic era: 100 MeV > T > me ( 200 ^ T > 1) 

The constituents of the cosmic fluid were leptons (neutrinos-antineutrinos 
and electrons-positrons), baryons (neutrons and protons), photons and cold 
dark matter. We will assume besides these particles a homogeneous (time 
dependent) magnetic field. At T ~ 1 MeV the neutrinos decouple, how¬ 
ever for temperature values such that T > me, {me — 0.5 MeV), pho¬ 
tons, neutrinos and electron-positrons continue with a unique temperature 


Te± =T^ = T, = T. 


— Epoch II: Beginning of the radiation epoch: me > T > 0.06 MeV (1 > 


r > 0.12) [SU-ISS] 


When the temperature drops below me, the electron-positron pairs are 
transformed into photons but not into decoupled neutrinos. After these 
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annihilations the number of photons is therefore greater than the number of 
neutrinos. Since thermal equilibrium is maintained until very few electrons 
remain, the entropy of the photon-electron-positron system for T > rrie 
is nearly the entropy of the photon system for T < mg. Thus, after the 
electron-positron annihilation the photon and neutrino temperatures are 
related by: 

T, = (4/ll)^/^T.,. (46) 

Since r ~ 0.029 s“^ x (t — ti), the cosmic times for these epochs range from an 
initial time at the outset of the leptonic era ~ 1 s = 0 at temperature 
T = 100 MeV, towards the final stage at the beginning of nucleosynthesis at 
t ~ 200 s ^ r ~ 6. 

Having defined the cosmic eras we are interested in studying, we undertake 


in this section the numerical integration of the system (28|-(32) and (33) 


under initial conditions specified at an initial time r = and following the 
assumptions summarised below: 

— An initial temperature of Tiji) ~ 200 and ending our analysis when tem¬ 
perature drops below T(ri) ~ 0.12. 

— Initial magnetic field is 0 ^ ^ (in Gauss 0 < B{ti) ~ 5 x 

lO^^G). 

— Initial shear is zero: 5'i(ri) = S 2 {Ti) = S' 3 (ri) = 0. From ([^-(30) and 


(37) this choice implies that all the anisotropy introduced by a Bianchi I 


geometry can be ascribed exclusively to the magnetic interaction through 
V± ^ 7^11. In other words, we have a magnetised Bianchi Universe that 
approaches an FLRW as the magnetic field becomes weak or negligible, 
while the strict limit B — 0 {B = 0) yields a pure FLRW evolution (not a 
perfect fluid Bianchi I model). While the shear tensor in early cosmic times 
should be absolutely negligible but not strictly (mathematically) zero, the 
assumption S'a(Ti) = 0 for a = 1, 2, 3 is a reasonable approximation. 

— The magnetic field modifies cosmic dynamics by both, the pure magnetic 
field contribution (Maxwell term) and its inclusion in the statistical treat¬ 
ment of the electron-positron gas. While the Maxwell term is present during 
the whole evolution, the electron-positron gas is considered as a magnetised 
gas only during the leptonic era. 


7.2 Kinematics 

We examine in figure the local proper volume for different initial values B(ri) 
of the magnetic field. The figure reveals a faster rate of expansion for larger 
B{Ti). However, if we consider values of B{Ti) compatible with observational 
bounds [28j . then the volume expansion is practically coincident with the ex¬ 
pansion rate of the FLRW model obtained when we set B{Ti) = 0 (see the 
grey line and the black triangles in the enclosed graph). In other words, the 
anisotropy effects on the volume expansion are completely negligible when the 
magnetic field complies with observational bounds. Since we are assuming zero 
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T = H.(t-t.) T = 


Fig. 1 Proper volume as function of time for increasing values of the initial magnetic field. 
The panel in the left hand side shows that larger values of Biji) produce faster expansion 
rates, while the curves of the panel in the right hand side show that the expansion rate of 
V(t) is practically indistinguishable from an FLRW expansion if we consider values of B{Ti) 
that fit observational constraints m- 


initial shear, the magnitude of the initial magnetic field intensity B{Ti) should 
determine the difference of the expansion rates for the three metric coefficients 
oi, 02 , 03 . This can be appreciated in figure which displays these functions 
for various values of B{Ti), showing distinct curves for these metric functions 
{i.e. larger anisotropy) for larger values of B{Ti) and an almost isotropic FLRW 
expansion (almost the same evolution for the three scale factors) for a weak 
field. 


7.3 Thermodynamics 

The anisotropy inherent in the magnetic field can also be appreciated through 
the anisotropy of the pressures. The classical Maxwellian contribution B^ pro¬ 
duces a decelerated fluid expansion in the direction of the magnetic field (z), 
as the latter produces a negative pressure or tension in its direction, and an 
accelerated expansion in the directions perpendicular to it due to a positive 
pressure. However, we obtain the opposite effect from the statistical contri¬ 
bution of B in the equation of state for the magnetised electron-positron gas 
mixture: the pressure in the direction parallel to the field increases and that 
in the perpendicular direction decreases [39ll^l46l[5T] . Although the dominant 
effect in the dynamics comes from the contribution B^, we can observe that 
including the magnetised electron-positron gas leads to a slight acceleration of 
the cosmic rate of expansion. In fact, the contribution of the electron-positron 
gas is of the order of 0 {a) x B"^ (where a is the constant of fine structure) [SS] , 
therefore this contribution will always be sub-dominant in comparison with 
that of the “pure” magnetic field ~ B^. 
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Fig. 2 Metric coefficients ai(r), a. 2 (T), a 3 (T) for varied values of Biji). Larger values of 
B{Ti) lead to very different curves (large anisotropy) for the three scale factors, while smaller 
values lead to almost the same curves for them (almost isotropy). 


The dominance of the Maxwellian contribution in the fluid expansion 
is consistent with the curves displayed in figure the scale factor 03 (r) (di¬ 
rection parallel to B) expands at a slower rate than 02 ( 1 ") and ai(T) (direc¬ 
tions perpendicular to B). This dominance is also displayed in figures and 
where we considered magnetic held values much larger than observational 
bounds to highlight this effect. Figure shows that the pressure V± is much 
smaller than V \\, while hgure depicts the growth of the volume expansion 
for different interactions of the magnetic held: the electron-positron as the 
magnetised/ ( non-magnetised ) gas (solid / ( dashed ) curves) and the grey curves 
(solid or dashed) represent the volume expansion obtained after eliminating 
the Maxwell term in the pressures and energies. 

So far, we have not mentioned the contribution of the magnetic held to the 
energy density and to the huid dynamics. This contribution depends on its ini¬ 
tial value B{Ti). It is negligible when B{Ti) T^{Ti) but becomes signihcant 
when B{Ti) ~ T^{Ti). In the former case B and T exhibit almost identical be¬ 
haviour to that of the well known FLRW solutions (along the lines of [S]) that 
follow by assuming T ~ 1/a and B ^ 1/a^, where a is the FLRW scale factor 
for a radiation dominated huid. However, if B{Ti) ^ T'^{Ti) there are signihcant 
diherences with the FLRW evolution. In order to explore these diherences we 
depict in Figure]^ the energy density for the different components of the huid 
mixture during the leptonic era, assuming an (unrealistic) initial magnetic held 
B{Ti) = 2 X T^{Ti) that is much larger than values allowed by observational 
bounds. As follows from these graphs, the curve that corresponds to the pure 
magnetic Maxwell term contribution (oc B^) decays very fast, but is overtaken 
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Fig. 3 Perpendicular and parallel pressures for each mixture component. We assumed 
B{Ti) = 2 X T^(ri). See further explanation in the text. 



T = 


Fig. 4 Contribution of each magnetic field effect in the volume expansion. We assumed 
B(r,) = 2x7^(ri). 


by the curves for the neutrinos, the electrons-positrons and the photons. This 
behaviour is very different from the expected radiation-like FLRW behaviour 
T ^ 1/a and B ^ 1/a^ of much weaker magnetic fields (subjected to ob¬ 
servational constraints). Hence, if B{Ti) ~ the Bianchi I and FLRW 

dynamics lead to very different decay rates for both the temperature and the 
magnetic field (see Sec. in particular equations ( |44| - @). 
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Fig. 5 Energy densities (total and for each mixture component). We assumed B{Ti) = 
2 X T\n). 


7.4 The magnetic field 


The different behaviour of the pure magnetic energy density (oc is il¬ 
lustrated in figure The solid curves represent the solutions obtained from 
Bianchi I dynamics, whereas the dashed curves correspond to a magnetic field 
in an FLRW context {B oc 1/a^). In both cases we plotted in the left hand side 
panel curves for two values of the initial magnetic field: B{Ti) = 27^(r^) (black 
line) and B{Ti) = T^(ri)/2 (grey line). In the right hand side panel we plot¬ 
ted the relative error Z\g(r) between B{t) obtained from Bianchi I dynamics 
and its FLRW equivalent B{Ti) / {t) for weak magnetic fields complying with 
{Bq) < 10“®G (solid line) and (i?o) I0“®G (dashed line). As shown by this 

graph, Ais(t) increases as the fluid evolves. This is an interesting result: while 
the evolution of kinematic and state variables for initial values Blji) of weak 
fields is practically indistinguishable in Bianchi I and FLRW dynamics, the 
evolution of the magnetic field itself can reveal non-negligible differences. The 
explanation for this effect is straightforward: while « 0102 holds for weak 
fields, each Bianchi I scale factor is related to the FLRW scale factor by a small 
(but time dependent) correction: aiir) ~ a(r)-|-ei(r) and 02 (r) ~ a(T)-|-e 2 (T), 
hence we have: 


[B] 


[^]flrw 


{€i{t) + e2{T))Bi 
a^{T) 


(47) 


where and [B]^^ are the FLRW and Bianchi I scaling laws given by ( |43 

and (44). Since we are assuming that the shear tensor vanishes at r = r^, then 
= e 2 (Ti) = 0, and thus the Bianchi I and FLRW forms initially coincide. 
However, as the expansion proceeds the “error” introduced by ei(r) -|- e 2 (T) 
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T=H.(t-t) 



T = H.(t-t) 


Fig. 6 Pure magnetic energy in the magnetised Bianchi I model and in a FLRW con¬ 
text (assuming B oc a“^) for various values of Biji) (left hand side panel). Relative error: 

= ^(t) — B{Ti)/a?[T)\/B{T) is plotted in the right hand side panel for magnetic 
fields complying (Bo) $ 10“®G (solid line) and (Bo) $ 10“®G (dashed line). 


in (47) (which enters in Z\b(t)) will be negligibly small only for 


but 


necessarily grows to small (but not necessarly negligible values) as the fluid 
expands. Although the relative error remains negligible for the calculations 
using the more restrictive bound {{Bq) < 10“®G), is of the order of 10“"^, 
it and can reach values of the order of 10“^, if we consider (Bg) < 10“®G. On 
the other hand, since even with these corrections the field itself is very weak, 
the evolution of the kinematic and state variables is practically insensitive to 
them. 


Finally, we remark that larger initial values of the magnetic field affect 
the dynamics of the other thermodynamical functions. Since the latter depend 
on the temperature, this can be illustrated by its different evolution with 
and without magnetic field. Hence, we depict in figure the evolution of the 
quantity 


Aj- = 


^LRW 

^LRW 



i^FLRW 


-T' 


T,- 


(48) 


which provides an estimation of the relative difference between the temper¬ 
ature in the magnetised Bianchi I mixture and the temperature in the non- 
magnetised FLRW model (7 (lrw) that results by setting B = Q. Since the 
initial temperature value is fixed, all curves depicted in the figure start at 
Aj- = 0. As T advances the evolution of the temperature of the magnetised 
mixture differs slightly (depending on the value of B{Ti)) from the FLRW evo¬ 
lution without magnetic field, with Aj- reaching a maximum value close to 
(when the magnetic field is stronger) and dropping as the expansion proceeds 
and the magnetic field decays. However, for weak helds values compatible with 
observational bounds Aj- ^ 0(10“®), therefore T is practically indistinguish¬ 
able from Tfi^rw 
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Fig. 7 The figure displays the relative ratio A defined in 
a magnetised Bianchi I fluid mixture and the same fluid 
FLRW metric. 


48| | between the temperature for 
with zero magnetic field and an 


8 Conclusions 

We have examined the evolution of a magnetised cosmic fluid mixture between 
the leptonic and nucleosynthesis eras by means of a non-perturbative approach 
to the Einstein-Maxwell field equations. For this study we considered a Bianchi 
I model, as it provides one of the simplest geometries that is compatible with 
the inherent anisotropy of the magnetic field. We considered as field sources a 
mixture of baryons (neutrons and protons), leptons (neutrinos, antineutrinos, 
electrons and positrons), cold dark matter WIMP’s and photons, together with 
an already existing time dependent magnetic field directed along the z-axis (we 
also assumed an nonzero cosmological constant). For the dynamical study of 
this cosmic fluid we transformed the Finstein-Maxwell equations into a first 
order system amenable for numerical integration. 

The general features of the behaviour of the physical quantities during 
the two epochs we are concerned can be summarised as follows. Temperature, 
magnetic field and thermodynamical functions depending on them, like pres¬ 
sures and energy densities of each species, start decaying at a fast rate and 
later at slower rates as the expansion proceeds. In epoch I at the end of the 
leptonic era the neutrinos provide the dominant contribution to the energy 
density and pressures, followed by the contribution of the electron-positron 
gas and photons. On the other hand, after the electron-positron annihilation 
the radiation dominates cosmic dynamics. These results agree with well known 
previous results [541155] . 

We have found that cosmic dynamics in the epochs described above is 
modified by both, the Maxwell term (~ B"^) and the presence of the magnetic 
field in the equations of state of a magnetised electron-positron gas. For weak 
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magnetic fields compatible with observational constraints these modifications 
are negligible (though this is not true for the evolution of the held itself). For 
magnetic helds complying with eB[ti) ~ the main contribution to the 

dynamics comes from the Maxwell term, leading to an anisotropic expansion 
that can be appreciated by the different time growth of the metric coefRcients 
and anisotropic pressures. The anisotropy caused by the Maxwell term pro¬ 
duces an accelerated expansion in the direction perpendicular to the held and a 
decelerated one in the held direction. On the other hand, the dependence on B 
in the electron-positron equation of state acts as a small correction 0 {a) x 
to the energy density where a is the hne structure constant (see further detail 
in [KH]b 

We found that the effect of weak magnetic helds compatible with obser¬ 
vational bounds on the evolution of the main thermodynamical variables are 
negligible. For such magnetic helds the time decay of these variables is prac¬ 
tically undistinguishable from their evolution under a pure FLRW dynamics. 
Since observational bounds yield a “near FLRW” Bianchi I model that is 
exactly spatially hat, it resembles some sort of “exact” perturbation whose 
evolution seems to be comparable with that of linear perturbations of mag¬ 
netised sources on a spatially hat FLRW background |23]. Lfowever, because 
of its homogeneity, the magnetic held in the “near FLRW” Bianchi I model 
permeates the whole time slices and thus cannot describe the scale depen¬ 
dence of standard cosmological perturbations, and thus cannot be used (even 
approximately) to examine early Universe magneto-genesis that involve seed 
helds undergoing a scale dependent amplihcation process (as for example the 
inverse cascade mechanism or cosmic plasma instabilities) to account for ob¬ 
served magnetic helds of astrophysical interest (see [Sl ITUllTnfHlfTSlfTB] 1. The 
only possibility (which was considered in [37]) is to regard the “near FLRW” 
magnetised Bianchi I model as providing an approximated description of mag¬ 
netic helds coherent in very large supra-horizon scales, which may be compat¬ 
ible with some inhationary magneto-genesis scenarios (see pununmi]). 

Given the limitations of the Bianchi I model mentioned above, we have 
provided a comparison of our model with the approach of isiiiiiiiHiiini, who 
considered a similar Stetistical Mechanics framework and examined magnetic 
helds during nucleosynthesis under a purely FLRW framework based on the 
assumption that the spacetime geometry is not affected by magnetic helds in 
these cosmic times (hence they rule out the anisotropy in the momentum- 
energy tensor that includes the magnetic interaction). While the results of jS] 
[371I331I33] can always be obtained from our approach under these same assump¬ 
tions (which may be well justihed for weak helds in a realistic cosmic huid), 
our non-perturbative methodology allows us to examine magnetised cosmic 
huids also under conditions in which such assumptions may not be justihed. 

A natural continuation of this study is revisiting the magnetic held con¬ 
straints as well as the study of abundance of light elements, all under the 
framework of Bianchi I dynamics applied specihcally to those magneto-genesis 
scenarios compatible with supra-horizon scale magnetic helds. This will be at¬ 
tempted in future work. 
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A Local kinematic variables 

The kinematics of local fluid elements can be described through covariant objects defined by 
the 4—velocity field . For a Kasner metric in the comoving frame endowed with a normal 
geodesic 4-velocity, the only non-vanishing kinematic parameters are the expansion scalar, 
0, and the shear tensor 


^ ^ ;CK 5 2 ^a:/3 ! 

where = UaU^ -\-9ap is the projection tensor and rounded brackets denote symmetriza- 
tion. For the Kasner metric these parameters take the form: 


ai a2 az 

cr^“ = diag [c7^, , ct/ , O] = diag [Ei , 1:2, ^3,0] , 

where: 

^ _ 2 ha _ 1 hb _ 1 

3 da 3 db 3 dc 

a^bT^c (a,b,c = 1,2,3). 


(50) 

(51) 


(52) 


The geometric interpretation of these parameters is straightforward: 9 represents the isotropic 
rate of change of the 3-volume of a fluid element, while cr^ describes its rate of local de¬ 
formation along different spatial directions given by its eigenvectors. Since the shear tensor 
is traceless: = 0, it is always possible to eliminate any one of the three quantities 

(Z*!, 172, ^ 3 ) in terms of the other two. We choose to eliminate Ui as a function of (172, 

Finally, the substitution of | |50[ l and into the Einstein equations, allow us to obtain 
a first order system of differential equations. 


B Particles 

At temperature values in the range 100 MeV > T > me neutrons and protons are free in 
chemical equilibrium. The equilibrium is possible because of reactions transforming neutrons 
into protons and vice versa m- 

UeUe~p, VeP<r^e'^n, (53) 

which implies the following variation of their numbers: 

Un _ 

— = e T , AM = mn — mp . (54) 

rip 
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On the other hand, for temperature values less than me the neutrons decay freely. Hence 
the densities of neutrons and protons satisfy 


hn = —Qrin — Prin . 


(55) 


hp = —Oup + Frin , 


(56) 


where F = llrn is the decay rate (rn 5 ^ 881.5 s). Also, from the charge neutrality we have 
n^-=np. 

Finally, since we assume during the whole evolution a baryons/photon ratio 77 = 5x 10“^, 
we have 


'^71 


(57) 


where rin, rip and are respectively the numbers density of neutrons, protons and photons. 
In this way we can roughly estimate the neutron and proton concentrations, a necessary task 
to obtain their rest energy. Since we are not interested in calculating element abundances, 
this rough approximated result is sufficient for our purposes. To improve these calculations 
the magnetic field should be included in the analysis [ 8 ], but this is beyond the scope of the 
present paper. 


B.l Nucleosynthesis constraints on light elements 


Since standard calculations of cosmological nucleosynthesis assume an FLRW universe with 
a radiation equation of state, it is worthwhile commenting on the effects of considering 
an anisotropic universe model on the primordial production of "^He. As discussed in [57], 
the time dependence of the radiation density is very important in determining the helium 
abundance. This yields the following bound for the shear eigenvalues in a Bianchi I model 
during the nucleosynthesis: Yp < 0.26 requires (cr/H) < 0.2 . In the previous equations 
Yp and H denote the primordial mass fraction of '^He and expansion scalar, respectively, 
and cr^ = + F 2 + ^ 3 ^) /2. However, for the above considered models with magnetic 

fields such that (Bq) < I0~®G, we have that {a/H) <C I remains valid throughout the 
nucleosynthesis process (notice that we assumed the anisotropy of the fluid to be caused 
only by the magnetic field). 


References 

1. R. Wielebinski, Magnetic Fields in the Milky Way, Derived from Radio Continuum Obser¬ 
vations and Faraday Rotation Studies. In: Wielebinski R, Beck R (eds) Cosmic Magnetic 
Fields, Lecture Notes in Physics, Berlin Springer Verlag, vol 664, p 89 (2005). 

2. P. P. Kronberg, Rept. Prog. Phys. 57, 325 (1994); R. Beck, Space Science Reviews 166, 
215 (2012). 

3. P.P. Kronberg, J.J. Perry, E.L.H. Zukowski Ap.J. 387, 528 (1992); M.L. Bernet, F. 
Miniati, S.J. Lilly, P.P. Kronberg, M. Dessauges-Zavadsky, Nature 454, 302 (2008); A. M. 
Wolfe, R.A. Jorgenson, T. Robishaw, C. Heiles, J.X. Prochaska, Nature 455, 638 (2008). 

4. T. E. Clarke, P. P. Kronberg and H. Boehringer, Astrophys. J. 547, LIII (2001) 
C. L. Carilli and C. B. Taylor, Ann. Rev. Astron. Astrophys. 40, 319 (2002) C. Vogt 
and T. A. Ensslin, Astron. Astrophys. 434, 67 (2005) A. Bonafede, L. Feretti, M. Murgia, 
F. Covoni, C. Ciovannini, D. Dallacasa, K. Dolag and C. B. Taylor, Astron. Astrophys. 
513, A30 (2010) 

5. E.J. Kim, P.P. Kronberg, C. Ciovannini and T. Venturi, Nature 341, 720-723 (1989); 
Y. Xu, P. P. Kronberg, S. Habib and Q. W. Dufton, Astrophys. J. 637, 19 (2006) 

6 . P. P. Kronberg, Rep. Prog. Phys. 57, 325 (1994); R. Beck, A. Brandenburg, D. Moss, 
A. Shukurov and D. Sokoloff, Ann. Rev. Astron. Astrophys. 34, 155 (1996); F. Covoni 
and L. Feretti, Int. J. Mod. Phys. D 13, 1549 (2004) 




26 


I. Delgado Caspar et al. 


7. A. Neronov and 1. Vovk, Science 328, 73 (2010) F. Tavecchio, G. Ghisellini, L. Foschini, 
G. Bonnoli, G. Ghirlanda and P. Coppi, Mon. Not. Roy. Astron. Soc. 406, L70 (2010) 
K. Dolag, M. Kachelriess, S. Ostapchenko and R. Tomas, Astrophys. J. 727, L4 (2011) 

8. D. Grasso and H. R. Rubinstein, Phys. Rept. 348, 163 (2001) 

9. A. Kandus, K. E. Knnze and C. G. Tsagas, Phys. Rept. 505 (2011) 1 

10. L. M. Widrow, D. Ryu, D. R. G. Schleicher, K. Subramanian, C. G. Tsagas and 
R. A. Treumann, Space Sci. Rev. 166, 37 (2012) 

11. R. Durrer and A. Neronov, Astron. Astrophys. Rev. 21, 62 (2013) 

12. J. Martin and J. Yokoyama, JCAP. 1, 25 (2008); K. Subramanian, Astron. Nachr. bf 
331, 110 (2010); C. Caprini and L. Sorbo, JCAP. 1410, 056 (2014) 

13. J.D. Barrow and C.G. Tsagas, Phys. Rev. D 77 107302 (2008) 

14. G. B. Field and S .M. Carroll, Phys. Rev. D 62, 103008 (2000); R. Banerjee and 
K. Jedamzik, Phys. Rev. D 70, 123003 (2004) 

15. P. Arnold and G. D. Moore, Phys. Rev. D 73, 025006 (2006); P. Arnold and G. D. Moore, 
Phys. Rev. D 73, 025013 (2006); P. Arnold, Int. J. Mod. Phys. E 16, 2555 (2007); P. Arnold 
and G. D. Moore, Phys. Rev. D 76, 045009, (2007); P. Arnold and P. S. Leang, Phys. Rev. 
D 76, 065012 (2007) 

16. K. Ichiki, K. Takahashi, H. Ohno, H. Hanayama, and N. Sugiyama, Science, 311, 827 
(2006); S. Saga, K. Ichiki, K. Takahashi and N. Sugiyama, Phys. Rev. D 91, 123510 (2015) 

17. A. Sanchez, A. Ayala and G. Piccinelli, Phys. Rev. D 75, 043004 (2007); 

18. J. Navarro, A. Sanchez, M. E. Tejeda-Yeomans, A. Ayala and G. Piccinelli, Phys. Rev. 
D 82, 123007 (2010); 

19. G. Piccinelli, A. Sanchez, A. Ayala and A. J. Mizher, Phys. Rev. D 90, no. 8, 083504 
(2014) 

20. C. G. Tsagas, A. Challinor and R. Maartens, Phys. Rept. 465 (2008) 61 

21. J. D. Barrow, R. Maartens and C. G. Tsagas, Phys. Rept. 449, 131 (2007) 

22. C. Tsagas and R. Maartens, Phys. Rev. D 61 (2000) 083519 

23. C. G. Tsagas and R. Maartens, Class. Quant. Grav. 17 (2000) 2215 

24. C. A. Clarkson, A. A. Coley, R. Maartens and C. G. Tsagas, Class. Quant. Grav. 20 
(2003) 1519 

25. C.G. Tsagas and J.D. Barrow, Class. Quantum Grav. 14, 2539 (1997); C.G. Tsagas 
and J.D. Barrow, Class. Quantum Grav. 15, 3523 (1998); C.G. Tsagas and R. Maartens, 
Phys. Rev. D 61, 083519 (2000); 

26. C. Zunckel, G. Betschart, P. K. S. Dunsby and M. Marklund, Phys. Rev. D, 73, 103509 
(2006); B. Mongwane, P. K. S. Dunsby and B. Osano, Phys. Rev. D 86, 083533 (2012) 

27. D. Grasso and H. R. Rubinstein, Astropart. Phys. 3, 95 (1995) 

28. D. Grasso and H. R. Rubinstein, Phys. Lett. B 379 (1996) 73 

29. M. Kawasaki and M. Kusakabe, Phys. Rev. D 86, 063003 (2012) 

30. H. J. Hortua, L. Castaneda, and J. M. Tejeiro, Phys. Rev. D, 87, 103531 (2013). 

31. Zeldovich, Ya.B., Zh. Eksp. Teor. Fiz. 48, 986 (1964) (English translation in JETP, 21, 
656 (1965)). 

32. K. S. Thorne, Astrophys. J. 148, 51 (1967). 

33. V.G. LeBlanc, Class. Quantum Grav. 15, 1607 (1998); V.G. LeBlanc, D. Kerr and 
J. Wainwright Class. Quantum Grav. 12, 513 (1995); J. T. Horwood and J. Wainwright, 
Gen. Rel. Grav. 36, 799 (2004); I. Singh Kohli and M. C. Haslam, Phys. Rev. D 88, 
063518 (2013) 

34. V.G. LeBlanc, Class. Quantum Grav. 14 2281 (1997) 

35. K. C. Jacobs, Astrophys. J. 155, 379 (1969). 

36. E. J. King and P. Coles, Class. Quant. Grav. 24, 2061 (2007) 

37. J. Adamek, R. Durrer, E. Fenu and M. Vonlanthen, JCAP 06 017 (2011) 

38. T. Harko, F. S. N. Lobo and M. K. Mak, Galaxies 2, 496 (2014); P. Fleury, C. Pitrou and 
J. P. Uzan, Phys. Rev. D 91, no. 4, 043511 (2015); T. Schucker, A. Tilquin and G. Valent, 
Mon. Not. Roy. Astron. Soc. 444, no. 3, 2820 (2014); P. A. R. Ade et al. [Planck Collab¬ 
oration], Astron. Astrophys. 571, Al (2014); P. A. R. Ade et al. [Planck Collaboration], 
Astron. Astrophys. 571, A15 (2014); P. A. R. Ade et al. [Planck Collaboration], Astron. 
Astrophys. 571, A16 (2014); E. Komatsu et al. [WMAP Collaboration], Astrophys. J. 
Suppl. 192, 18 (2011) 

39. M. Chaichian, S. S. Masood, C. Montonen, A. Perez Martinez and H. Perez Rojas, Phys. 
Rev. Lett. 84, 5261 (2000) 



A non-perturbative study of the evolution of cosmic magnetised sources 


27 


40. A. P. Martinez, H. P. Rojas and H. J. Mosquera Cuesta, Eur. Phys. J. C 29, 111 (2003) 

41. E. J. Ferrer, V. de la Incera, J. P. Keith et al.^ Phys. Rev. C 82, 065802 (2010) 

42. Landau, L. D., and E. M. Lifshitz. The Classical Theory of Fields Addison-Wesley. 
Reading Mass, and Pergamon, London (1971). 

43. A. Ulacia Rey, A. Perez Martinez and R. A. Sussman, Gen. Rel. Grav. 40, 1499 (2008) 

44. A. Ulacia Rey, A. Perez Martinez and R. A. Sussman, Int. J. Mod. Phys. D 16, 481 
(2007) 

45. D. Manreza Paret, A. Perez Martinez, A. Ulacia Rey and R. A. Sussman, JGAP 1003, 
017 (2010) 

46. 1. Delgado Gaspar, A. Perez Martinez, R. A. Sussman and A. Ulacia Rey, European 
Journal of Physics G, 73, 2502, (2013) 

47. J. D. Barrow, Phys. Rev. D 55, 7451 (1997) 

48. J. D. Barrow, P. G. Ferreira and J. Silk, Phys. Rev. Lett. 78, 3610 (1997) 

49. E. N. Parker, Cosmological Magnetic Field, Oxford, Clarendon (1979). 

50. L. Mestel, Stellar Magnetism, Oxford, Oxford University Press (1999). 

51. H. Y. Chiu, V. Canute, and L. Fassio-Canuto, Phys. Rev. 176 (1968) 1438; V. Canute 
and H. Y. Chiu,Phys. Rev. , 173,1229, (1968); V. Canuto and H. Y. Chiu, Phys. Rev. , 
173, 1220, (1968); V. Canuto, and H. Y. Chiu, Phys. Rev. , 173, 1210, (1968). 

52. J. Schwinger, Phys. Rev. 82, 664, (1951) . 

53. G. Greenstein, Nature, 223, 938, (1969). 

54. S. Weinberg, Cosmology, Oxford University Press Inc, New York (2008). 

55. Rich, James. Fundamentals of cosmology. Springer Science Business Media, (2009). 

56. P. J. Kernan, G. D. Starkman and T. Vachaspati, Phys. Rev. D 54, 7207 (1996) 

57. Wainwright, John, and George Francis Rayner Ellis. Dynamical systems in cosmology. 
Cambridge University Press, (2005). 



